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r 3 1 | Abstract. We consider quantum Hamiltonian systems composed of mutually 

q ■ interacting "dynamical subsystem" with one or several degrees of freedom and 

'thermostat" with arbitrary many degrees of freedom, under assumptions that the 
£^ interaction ensures irreversible behavior of the dynamical subsystem, that is hnitc 

diffusivities of its coordinates in thermodynamically equilibrium state and finite drift 
velocities and mobilities in non-equilibrium steady state in presence of external driving 
^ " forces. It is shown that, nevertheless, regardless of characteristics of the interaction, the 

diffusivity and mobility have no certain values but instead vary from one observation 
to another and undergo 1/f- type or flicker- type low-frequency fluctuations. 

a 

O . PACS numbers: 05.30.-d, 05.40.-a, 05.60.Gg, 71.38.-k 

Keywords : Dynamical foundations of kinetics, Quantum Brownian motion, Quantum 
transport, Quantum kinetic equations, Mobility fluctuations, 1/f noise, Low-frequency 
^ \ noise, Flicker noise 

q 

^! 1. Introduction 

(N 

Nearly 30 years ago in [1, 2] and then in [3, 4, 5, 6] an original explanation of electronic 
1/f-noise was suggested by Prof. G. N. Bochkov and me and putted in phenomenological 
^ ■ statistical model which easy produces correct estimates of the noise level and directly 

connects to rigorous statistical mechanics (for later looks at this model see [7, 8]). Key 
idea of our explanation was that own (thermal) fluctuations in rates of various transport 
processes in many-particle systems do not cause a "back reaction" and therefore have no 
definite "relaxation times" . Hence, that are scaleless fluctuations, with power-law low- 
frequency spectrum (which thus has no relation to folklore "composition of lorentzians" ) . 

In [1, 2, 3, 4, 5, 6] this principal idea was exploited in specific terms of electric charge 
transport and "Brownian motion" of charge carriers (one may see also [9, 10, 11, 12]). 
But, clearly, it can be applied also to to different processes in large variety of physical 
systems (in most wide sense of the word "physical" [7]). Particularly, to Brownian 
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motion ("random walks") of a fluid particles [13, 14, 15, 16, 17, 18, 20], as well as to 
collective fluctuations and non-equilibrium processes in fluids [15, 16] and solids [21]. 

Simultaneously, the works [13, 21] and later [14, 15, 16, 17, 18] and [11], in fact, 
presented substantiations of the mentioned key idea on the base of rigorous statistical 
mechanics, and confirmed the sentence from [5] that "1/f- noise is a kind of tribute to be 
paid to the dynamics for dissipation and irreversibility properties of physical systems" 
(my translation from Russian original). 

This sentence fully agrees with results of the deep N. Krylow's investigation [22] 
(discovered by me when preparing [13]) where it was shown that in many-particle 
statistical mechanics, firstly, relative frequencies of observable events can be different at 
different system's phase trajectories (experiments), regardless of duration of observations 
and time averaging (see characteristic mental example in introductions in [16] and [14] 
and more examples in [7]). Therefore, generally one has no rights to presume any a 
priori certain (even let unknown numerically) "probabilities" of events. Secondly (as 
the consequence), it is possible that statistical correlations do exist even betwen events 
what are independent in physical (cause-and-consequence) sense. 

In [13, 21, 14, 15, 16] I demonstrated that real "fundamental" 1/f-noise (i.e. such 
1/f-noise whose spectrum has no saturation at zero frequency) is just manifestation of 
the Krylov's uncertainty of events' probabilities (additional explanations can be found 
in [7] and introductory or discussion-resume sections in [8, 9, 17, 18, 19] and [11]). 

Unfortunately, the enumerated works, - including the Krylov's book, - yet has not 
excited any scientific-community's response except neglect or aversion f . Interestingly, 
such aversion was pointed out and commented already by Krylov himself [22] , as product 
of common prejudices educated by standard "probability-theoretical" view at physical 
world and its statistical description § . Hence, we have to continue activity in this field 
as far as possible. 

In this paper our aim is to illustrate the aforesaid ideology (and once again confirm 
its correctness) by example of simple Hamiltonian microscopic model of quantum 
Brownian motion. More concretely, to show that if an interaction of (quantum) 
particle with (quantum) thermostat make this particle "Brownian", that is ensures 
diffusive character of its chaotic motion, - when mean square of its displacement grows 

i 

Thus showing that progress in fundamental science also undergoes fundamental 1/f-type fluctuations. 
For example, manuscripts [15] were rejected by JSTAT and JSP, respectively, [16] by Physica A, [17] 
by JSP, EPL and JMP, [8] and [10]- [11] by JETP, [19] by PRE, JPA and JSP, etc., in all cases without 
meaningful or even any motivation. 

§ 

The modern notion of "dynamical chaos" still has not changed this situation, since, instead of 
investigating specificity of the chaos in (infinitely-) many-particle Hamiltonian systems, - where it 
naturally creates 1/f-noise [16, 21] and thus "kinetic non-ergodicity" [7, 16], - scientists try to squeeze it 
in "Procrustean bed" formed with "Marcovian partitions" , "Bernoullian flows" and primitive stochastic 
processes with a priori introduced "probabilities of elementary events" . Although there is also concept 
of such "K-flows" what can not be re-coded to Bernoullian flows [23]. 
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proportionally to time, - then such interaction ensures also non-zero 1/f- type fluctuations 
of its diffusivity (and hence mobility). In other words, the latter has no certain value 
but varies from one measurement to another. 



2. A class of systems under consideration 

Let us consider a quantum system with Hamiltonian 

H = H d (q,p) + H b (q) , (1) 

2 

H d (q,p) = ^- -f-q , (2) 
2m 

where Hd(q,p) represents "dynamical subsystem" with canonical variables (coordinates 
and momenta) q and p, whose operators satisfy standard commutation relations 

[Qa,Pfi] = ih$ap , (3) 

and H b (q) represents "thermal bath" ("thermostat") along with interaction between it 
and the dynamical subsystem (DS). 

Our first principal assumption is that our system is "translationally invariant" in 
respect to the coordinates q, in the sense that properties of the operator H b (q), - as 
considered in the thermostat's Hilbert (phase) space, - do not depend on parameters 
q. Formally, this means that there exist such Hermitian operators II, - defined in the 
thermostat's space, - that 

H b (q + a) = exp (— iall/h) H b (q) exp (iall/h) (4) 

In other words, from thermodynamical point of view the thermostat is indifferent to g's 
values. 

Under such assumption, one can treat variables (observables) q and v = p/m 
like coordinate and velocity of unrestricted (generally multi-dimensional) "Brownian 
motion" . 

We want to consider possible statistical characteristics of these motion, at that 
basing, of course, on the von Neumann (quantum Liouville) evolution equation for full 
system's density matrix p, that is 

p = [H,p]/tH = Cp , (5) 

with [A, B] = AB — BA and C being the Liouville super-operator. 



3. Quasi-classic representation 

It is convenient to consider Eq.5 in the form most closely unifying quantum and 
classical treatments of DS's variables, that is in quasi-classical or Wigner representation. 
Simultaneously, by reasons what will be clear later, we want to go from the density 
matrix, p, to an equivalent characteristic function of variables of our system. We will 
make this in three steps. 
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1. Let | q) be eigen- vectors of operators q, and let us introduce functions-operators 

p(t,x,y) = (x + y/2\p\x-y/2) , (6) 
p(t,x,p) = J exp(-ipy/h)p(t,x,y)d d y/(27ih) d , (7) 

where d is number of pairs {q,p}, i.e. degrees of freedom of the dynamical subsystem 
(DS). The function-operator p(t,x,p) is jointly the Wigner's probability distribution 
function of DS's variables, - with x representing the coordinates, - and density matrix 
of the thermostat. 

From viewpoint of this distribution, evidently, 

P(t,x,y) = J exp(ipy/h) p(t,x,p)d d p (8) 

is nothing but characteristic function of the momenta p, in the same sense of the words 
"characteristic function" (CF) as in the probability theory. Or, to be more precise, 
p(t, x, y) is a hybrid of CF, - in respect to DS's momenta, - and probability distribution 
(PD), or density matrix, in respect to DS's coordinates and all thermostat's variables. 

2. Next, let us make transform 

p(t,x,y) = exp(—ixU/h)R(t,x,y)exp(ixU/h) , (9) 
p(t,x,p) = exp (— ixll/h) R(t, x,p) exp (ixll/h) 

After it, in terms of new (distribution) function-operator (density matrix) R(t,x,y), 
Eq.5 takes form 

R=--f^- + C f R = CR, (10) 
m oxoy 

C f R= l -f-yR+-^-[U,R}+^-{ H b (y/2) R - RH b (-y/2) } (11) 
n m ay in 

Clearly, one can say that R describes thermostat in movable frame connected to DS 
coordinates q (via relative coordinates, if say in classical language). 

One more Fourier transform, in addition to (8), in respect to x, produces function- 
operator 

R(t,k,y) = J exp(ikx) R(t,x,y)d d x , (12) 

which is joint semiclassical CF of all DS's variables "hybridized" with thermostat's 
density matrix. The corresponding evolution equation directly follows from Eq.10 : 

• hk dR „ „ „ „ . . 

R= — —+C f R = CR 13 

m ay 

3. For last step, we assume that, - as usually in many-particle physics (see e.g. 
[24]), - all possible states of the thermostat can be described as results of actions of 
definite particles (or quasi-particles) and/or quanta creation and annihilation operators, 
c\ and c s , onto "vacuum state". Index "s" here enumerates various sorts and modes 
of particles and quanta. Then all thermostat's observables, including H b (q) and II, 
can be treated as functions of c\ and c s . Correspondingly, thermostat's component of 
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full density matrix of our system can be completely described in terms of quantum CF 
defined e.g. like Tr# exp (z* c s ) p exp (z s c\) , where Tr# denotes trace over thermostat 
space, the repeated indices "s" mean summation (or/and integration), and z s and z* are 
test, or probe, parameters, either complex c- numbers or Grassmann numbers, depending 
on whether c\, c s are Bose or Fermi operators (at that, generally, z s and z* are thought 
as independent, not mutually conjugated, variables). 

Then, analogously, let us introduce full CF of our system by 

F{t, k, y, z,z*} = Ty b exp (z* c s ) R(t, x, y) exp (z s c\) (14) 

Importantly, under the coherent-state representation of R such defined CF coincides 
with mere CF of corresponding "quasi-probability density". Therefore ^{t, k,y, z,z*} 
is fully semi-classical object, and we can speak about it applying usual "probability- 
theoretical" terminology. 

In order to write out an evolution equation for it, first for any operator A, - 
composed of various c\ and c s , - let us introduce two operators A + and A~ acting 
in space of functions of the probe parameters and defined as follow: 

exp (z s 4) exp (z* c s ) A = 
( d d \ 

= A + [z, z*, — , — exp (z s 4) exp (z* s c s ) , (15) 



dz ' dz* r 
A exp {z s c\) exp (z* s c s ) = 

= A' (^z, z*, ^, -^j exp (z s 4) exp (z* c s ) (16) 

Thus A + and A~ are composed of multiplications and differentiations in respect to 
probe parameters. Obviously, structure of A + and A~ is unambiguously determined by 
that of A and commutation rules of 's and c's. At that, anyway 

Exploiting these definitions, one easy transforms Eq.13 into equation for the full 
CF F{t,h,y, z,z*}: 



T = 



hk d 

— IT + C f 

m ay 



7 = CT , (18) 



Here V,^(q) and V ± are the operators produced, - according to Eqs. 15-16 , - from H b (q) 
and II, respectively, and their arguments are omitted for brevity. 

Further, let us discuss construction and properties of solutions of Eq.18 and 
equivalent Eqs. 10 and 13. 



4. Equilibrium and steady non-equilibrium states 



First, consider stationary solutions of Eq.18, 13, 10 and 5. Because of the thermostat's 
property (4) we can expect that stationary solutions of Eq.10 are independent on x , 
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i.e. "spatially uniform", if one interprets the DS as "Brownian particle" (BP), while 
the thermostat as a medium where BP moves (then property (4) means thermostat's 
spatial infiniteness and uniformity). Correspondingly, Eqs.13 and 18 have (non-trivial) 
stationary solutions only when the parameter k ("wave vector") equals to zero, k — 0. 
Denoting them by R s t(y] f) and T st = J- St {y, z, z*; /}, we can write 

C f R st = , LjTst = , (20) 

with Cf defined by (11) and (19), respectively. 

Among all possible solutions of this equation, most physically interesting are that 
characterized by definite thermostat's temperature. In absence of the "external driving 
force" , i.e. at / = , such solution is thermodynamically equilibrium one determined 
by 

ReM = PeM oc {y/2\ exp {-[p 2 /2m + H b (q)]/T} | - y/2) , (21) 
Feq{y, z,z*} = Tr B exp (z* c s ) R eq (y) exp (z s c\) , (22) 

along with obvious normalization condition Tr# R eq (0) = 1 , J- eq {0, 0, 0} = 1 . 

Further, let us take these equilibrium expressions to be initial conditions to Eqs.13 
and 18 with / = const ^ and k — . Then stationary asymptotic of their solutions 
at t — > oo will give us such non-equilibrium solution of Eq.20 which can be treated as 
perturbation of R eq and T eq and therefore also is characterized by definite thermostat 
temperature. 

5. Characteristic functional, statistical correlations, cumulants, and 
statistics of the Brownian motion 

Of course, we just have made second principal assumption, namely, that construction 
of the Hamiltonian H b (q) ensures existence of the mentioned stationary asymptotic, at 
least at sufficiently small finite |/| . 

If it is so, then we can investigate stationary random walk of the BP, - i.e. random 
changes of the DS variable q(t) , - by considering non- stationary solutions of Eqs.10 or 
13 or 18 at t > with initial conditions 

R(t = 0, x, y) = 8{x) R st (y; f) , R(t = 0, k, y) = R st (y; f) , 

F{t = 0, k, y, z, z*} = F st {y, z, z*; /} , (23) 

respectively. Evidently, these initial conditions are quasi-classical representations of such 
state of DS (BP) which is balanced in respect to its interaction with thermostat but, 
at the same time, localized in the q -space at q(t = 0) = 0. Therefore solution of Eq.10 
for t > determines probability density distribution, W(t,x), of x(t) = q(t) — q(0) 
(BP's path, or displacement, during time t), while Eqs.13 and 18 characteristic function 
of this distribution: 

W(t, x) = Ti B R(t, x, y = 0) = (S(x{t) - x)) , (24) 
W(t,k) = [ e lkx W(t,x)d d x = 
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= Tt b R(t, k, y = 0) = F{t, k, y = 0, z = 0, z* = 0} = ( e ikx{t) ) (25) 

We introduced the angle brackets as standard comfortable designation of statistical 
ensemble averaging. With this designation, in the quasi-classical language, we can write 

F{t,h,y,z,z*} = (exp[ikx(t)+i^v(t) + z s c* s (t) + z* s c s (t)}) , (26) 

where v(t) = p(t)/m = dx(t)/dt represents velocity of q(t) 's changes, so that 

x(t) = fv{r)dr , (27) 
Jo 

and £ = my/h. 

Expression (26) visually shows that F{t, k,y, z, z*} is full characteristic function 
of all variables of the whole system "DS plus thermostat" (or, to be be precise, 
characteristic functional, since for infinitely large thermostat c s and c* form continuum 
set of variables). Taking in mind general properties of characteristic functions in the 
probability theory [25], we can write also 

F{t, k, y, z, z*} = Wit, k) Q{t, k, y, z, z*} F st {y, z, z*; /} , (28) 

where the middle multiplier on the right, 6 , contains all cross-correlations between the 
path x(t) , on one hand, and v(t) , c s (t) and c*(t) , on the other hand. This means 
that 

e{t, k = 0,y, z, z*} = G{t, k, y = 0, z = 0, z* = 0} = 1 (29) 

In terms of cumulants, i.e. irreducible correlations, - to be designated by double angle 
brackets, - 

In 6{f,*,V, *,**}= E ^|r« a;B (*)»7 , (0» = 

n,l = l 

= E ^TTT f ■■■ f {{<r l )...v{r n )r 1 \t)))dr l ...dr n , (30) 

n,l = l J ° J ° 

lnW(t,k) = E^«^(*)» , (31) 

n = l 

where r](t) = i£v(t) + z s c*(t) + z*c s (t) . 

Before discussion of possible time behavior of the x(t) 's cumulants in (31) and 
thus statistics of the Brownian motion, we have to realize most principal property of 
the evolution operator C in Eq.18 determining the x(t) 's cumulants. Namely, to take 
into account that 

6. Spectrum of the Liouville super-operator is purely imaginary 

Let \I/ Q is complete set of mutually orthogonal eigen-vectors of the full system's 
Hamiltonian (1): H\ty a ) = E a \^f a ) . Then p a p = \ty a ) (ty p\ is complete set of 
orthogonal eigen-vectors of the Liouville, or von Neumann, super-operator £ in (5): 
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Cp a p = i((Ep — E a )/K)p a p . The orthogonality is understood in the sense of the natural 
"scalar product" (A, B) = Tr A^B = Tr D Tr B A^B . 

Hence, £ has purely imaginary spectrum. This statement equally comprises, of 
course, the evolution operators in Eqs.10, 13 and 18 which are equivalent representations 
of original Liouville operator. The same can be said about the operators £/ . 

In other words, these operators can not have real eigen- values or complex ones with 
nonzero real parts. This trivial truth of statistical mechanics is sufficient ground for 
fundamental theorem to be claimed in next section. 

Before it, notice, first, that in the representation defined by Eqs.9 and 10 the 
mentioned scalar product reads 



In the representations producing Eqs.13 and 18, evidently, the variable ("wave vector") 
k in fact plays role of passive parameter of the Liouville super-operator C . Therefore 
its eigen- vectors naturally divide into layers (subsets) which correspond to different k 's 
and can be marked by k (such treatment of C 's eigen- values will appear below). At 
that, the scalar product (32) reduces, - in case of Eq.13, - to 



in each separate layer (an equivalent formula for the case of Eq.18 looks rather 
cumbersome and will be written out elsewhere). 

7. Theorem about uncertainty of relaxation, diffusion and dissipation rates 

The general hope of conventional "kinetic theory", or "physical kinetics", etc., is that a 
proper interaction of DS (a part of much greater system) with thermostat (the rest of the 
system) can impel the first of them to irreversible and stochastic behavior characterized 
by well definite "kinetic coefficients", relaxation and dissipation rates, etc., at least under 
the "thermodynamic limit" (for infinitely large thermostat). In respect to presently 
considered class of systems, this assumption means that all the cumulant functions in 
(30) and (31) are fast enough decaying (integrable) functions of time differences t — Tj 
and Tj — Tfc , respectively. That is all terms of series (30) tend with time to finite limits, 
so that 




(32) 




(33) 



Q{t, k, y, z, z*} — > 0{oo, k, y, z, z*} ^ oo , , 
while all terms of series (31) are asymptotically linear time functions: 
In W(t, k) ->■ X(ik)t + const , 



(34) 



(35) 




(36) 



71=1 



A n (/) = lim«x B (f)»/f , 



(37) 
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where (v) = (dx(t)/dt) = (x(t)/t) is mean velocity vector and D is diffusivity tensor. 
For "good" enough thermostat and small force / one expects also that (v) = /if with 
\x being mobility tensor. 

In terms of the probability theory [25], Eq.35 says that the Brownian path x(t) 
behaves as a random process with independent increments and infinitely divisible 
probability distribution, and Eq.36 specifies that this is diffusive process with 
asymptotically Gaussian probability distribution. 

But, fortunately, these conventional assumptions can not be true. Indeed, if they 
were true then expressions (34)- (35) as combined with Eqs.18 and 28 would imply that 
asymptotically 



which would means that \(ik) is eigen-value of the Liouville super-operator. But it 
is certainly impossible since X(ik) by its definition inevitably has non-zero (negative) 
real component, ?R,\(ik) < at k ^ (while, to be recalled, £ 's spectrum is purely 
imaginary!). 

Consequently, the hypothetical linear asymptotic (35)-(37) is wrong, and in fact 



for some n > 1 (or n > 3 at / = 0), that is all particular increments of the path 
x(t) are essentially statistically dependent one on another (correlated with each other) 
regardless of time distance between them. 

At that, values n — 1 , or n < 3 at / = , are excluded from candidates to the 
super- linearity (as well as to sub-linearity) by our assumptions about existence of steady 
non-equilibrium state (described byj- st ) with finite mean ("drift") velocity (v) , which 
means finiteness of BP's diffusivity in equilibrium state (described by T eq ) at / = 
(this follows from the Einstein relation D = T\x easy provable for our systems too). 

Therefore, the crash of the "independence of increments" impliess that generally the 
super-linear cumulants' asymptotic (39) takes place at any n > 1 (or n > 3 , if / = ). 
Then, it remains to realize that this statement is equivalent to statement that mobility 
and diffusivity possess low- frequency fluctuations like "flicker noise" or 1/f -noise. 

In other words, the mobility and diffusivity (as well as related "kinetic 
characteristics" of x(t) 's interaction with thermostat) do not have certain values but 
change from one experiment (measurement) to another, with a spread what anyway 
hugely exceeds limits suggested by the "law of large numbers" . 

8. Conclusion 

We have expounded rather general theorem (or, strictly speaking, a "storage" of future 
rigorous theorem) stating that interaction of one Hamiltonian (sub-) system with 
another, - let serving as an arbitrarily large thermostat, - never can ensure well certain 



X(ik) 0{oo, k, y, z, z*} F st {y, z, z*; /} 
= £ 6{oo, k, y, z, z*} ^{y, z, z*; /} 



(38) 



«*»(*)» 



— > oo 



(39) 



t 
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quantitative characteristics of irreversible and dissipative behavior of the first of them. 
Instead, these characteristics undergo significant 1/f -type ("flicker" type) low-frequency 
fluctuations. Citing once again myself from [5], "1/f-noise is a kind of tribute to be paid 
to the dynamics for dissipation and irreversibility properties of physical systems" . 

Particular variant of this theorem was obtained in [11]. For better understanding 
its pre-history, physical meaning and possible applications, see also works referred there 
and in the Introduction above. Detailed elaboration and concrete applications of our 
present result will be considered separately. 
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